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Abstract 

Charged Higgs boson production in association with a top quark could be the first in- 
dication of the existence of Higgs particles. The Tevatron Run-II started data-taking 
in April 2001 at ^/s = 1960 GeV and could probe the existence of a charged Higgs 
boson beyond the current mass limit. We study the pp thU.^ production process 
with Monte Carlo simulations in HERWIG and PYTHIA, comparing expected cross 
sections and basic selection variables. 
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Abstract 

Charged Higgs boson production in association with a top quark could be 
the first indication of the existence of Higgs particles. The Tevatron Run-II 
started data-taking in April 2001 at y/s = 1960 GeV and could probe the ex- 
istence of a charged Higgs boson beyond the current mass limit. We study the 
pp tbH^ production process with Monte Carlo simulations in HERWIG 
and PYTHIA, comparing expected cross sections and basic selection variables. 

1. INTRODUCTION 

Charged Higgs bosons are predicted by non-standard models, for example Two-Higgs Doublet Models 
such as the Minimal Supersymmetric Standard Model (MSSM). Thus, their detection and the measure- 
ment of their properties (such as the mass which is not predicted by any model) play an important role 
in the investigation of an extended Higgs sector and in the understanding of the generation of particle 
masses. The current limit on the charged Higgs boson mass is set by the LEP experiments at 78.6 GeV, 
independent of the Higgs boson decay branching fractions (3. At the Tevatron, charged Higgs bosons 
could be discovered for masses well beyond this limit. 

If the charged Higgs boson mass mH± satisfies mH± < mt — nih, where mt is the top quark 
mass and rrib the bottom quark mass, it could be produced in the decay of the top quark t bH+. 
This so-called on-shell top approximation (qq, gg —>■ tt with t — > bH"*") was previously used in the 
event generators. Throughout this paper this process is denoted by pp — > tt — > tbH^. Owing to the 
large top decay width (Ft ~ 1.5 GeV) and because of the additional diagrams which do not proceed 
via direct tt production |2l|3l|, charged Higgs bosons could also be produced beyond the kinematic 
top decay threshold. The importance of these effects in the threshold region was emphasized in the 
previous Les Houches proceedings [T| and the calculations \2. 3 1 are implemented in HERWIG |5l |6| 
and PYTHIA 1 8 1 ^ . The full process is referred to as pp tbH^ . Examples of the graphs contributing 
to the pp tbH"*" process are IQJ : 

in /b 

b b \y 




The t-channel graph is one example of a diagram which does not proceed via tt production. This graph 
contributes to enhanced particle production in the forward detector region. 

A charged Higgs boson with mH± < TUt decays predominantly into a r lepton and a neutrino. 
For large values of tan (3 (> 5), the ratio of the vacuum expectation values of the two Higgs doublets, 
this branching ratio is about 100%. The associated top quark decays predominantly into a W boson or a 
second charged Higgs boson, and a b-quark. The reaction 

pp ^ tbH± (t ^ bW^) (H± ^ T^Ur) (2) 

is a promising channel to search for the charged Higgs boson at the Tevatron. Simulations are performed 
at the centre-of-mass energy y/s = 1960 GeV and for tan (3 = 20. 
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2. COMPARISON OF PRODUCTION CROSS SECTIONS 

The expected production cross sections are determined using HERWIG and PYTHIA simulations, and 
are shown in Fig. ^ The default mass and coupling parameters of HERWIG version 6.5 and PYTHIA 
version 6.2 are used. The cross sections depend strongly on the top decay width over the investigated 
mass range. For the top width, the Standard Model (SM) value Ft = 1.53 GeV is used at = 
210 GeV and the width is increased as a function of the charged Higgs boson mass to Ft = 1.74 GeV 
at mH± = 70 GeV in both generators. The production cross section in HERWIG is about a factor 2 
larger compared to PYTHIA which can be attributed to the default choices of the standard parameters. 
It is mostly driven by the different choice of the heavy quark masses entering the Higgs-quark Yukawa 
coupling. In PYTHIA a running b-mass is used at the tbH"*" -vertex. For mH± = 150 GeV the b-quark 
mass of 4.80 GeV is reduced to rrib = 3.33 GeV, while HERWIG uses mb = 4.95 GeV both in the 
kinematics and at the vertex. Other relevant parameters are the default Parton Distribution Functions 
(PDFs) and the coupling constants a and ag, as well as the scales used for evaluating the PDFs and 
couplings, which are not the same in the default setups of the two simulation packages. 

Tests comparing the total cross sections from HERWIG and PYTHIA for identical choices of all 
above parameters were performed and confirmed that the two implementations of the hard scattering 
matrix elements coincide numerically. In this study, however, we maintain the default configurations 
of the two simulation packages. Hence, differences in the various distributions may be taken as an 
indication of the theoretical systematic errors affecting this process. 

Fig. 1: Charged Higgs boson production cross section at = 
1960 GeV for tan/3 = 20. For the pp -> tbH* process (thick 
lines), the HERWIG expectation is larger by about a factor 2 com- 
pared to PYTHIA because of the different default setups as described 
in the text. The differences between the two PYTHIA curves for the 
pp tbH^ and pp tt ^ tbH* processes instead result from 
top decay width effects and because of the additional diagrams which 
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3. COMPARISON OF BASIC SELECTION VARIABLES 

At the parton level, several distributions of variables related to the event topology are compared between 
HERWIG and PYTHIA simulations. In addition, differences in the distributions between the pp 
ihBj^ process and the pp — > tt ^ tbH='= subprocess are demonstrated. Each comparison is based on two 
samples of 10,000 generated events. Effects of the different event fragmentation schemes in HERWIG 
and PYTHIA could influence the comparison and they are not considered here. The detector simulation 
of the Tevatron experiments would reduce further the sensitivity of these comparisons. Figure |2l shows 
the transverse momentum and pseudorapidity rj of the following particles: 

a), b) The b-quark produced in association with the in the pp — > tbH^ (dots) and pp — > tt — > tbH^*^ 

(solid line) processes in PYTHIA for mH± = 165 GeV. 
c), d) The b-quark produced in association with the H='= in the pp tbH^ process in HERWIG (dots) 

and PYTHIA (solid line) for mH± = 150 GeV. 
e), f) The b-quark from the top quark decay (t bW^) in the pp — > tbH^ process in HERWIG (dots) 

and PYTHIA (solid line) for m^± = 150 GeV. 
g), h) The r lepton from the decay in the pp — > tbH^ process in HERWIG (dots) and PYTHIA 

(solid line) for rn,jj± = 150 GeV. 
The differences in the pT and -q distributions are clearly visible in Figs.|2la) and b) between the processes 
pp tbH^ and pp — > tt ^ tbH^. The HERWIG and PYTHIA simulations show good agreement 
in the kinematic distributions of Figs. |2lc) to h) for both b-quarks and the r lepton. The decay of the r 
lepton is not considered here, but it should be noted that spin correlations must be taken into account in 
the study of the final state particles 19. 101 . 
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Fig. 2: Distributions of cliarged Higgs boson selection variables at the parton level for ^/s = 1960 GeV and tan /9 = 20 for 
HERWIG and PYTHIA. The variables are described in the text. In a) and b) the differences are mainly from top off-shellness 
effects. In c) to h) each pair of curves is normalised to an equal area. The error bars on the dots indicate the statistical uncertainty. 



4. CONCLUSIONS 

At Tevatron Run-II, about 1000 pp thU^ events per 1 fb~^ at ^/s = 1960 GeV could be produced 
for mH± = 100 GeV and tan (3 = 20, while about 100 events are expected for m^± = 150 GeV 
These expected event rates will strongly be reduced when selection criteria are applied to separate signal 
and background events. For the default choices of mass and coupling parameters in HERWIG and 
PYTHIA we observe significant differences in the simulated total cross sections. We have also studied 
the shape of basic selection variable distributions and found good agreement between the HERWIG and 
PYTHIA parton level predictions in the pp — > tbH^ process. In comparison with the pp — > tt — > tbH^ 
subprocess, which was used in previous HERWIG and PYTHIA versions, for mjj± > 160 GeV the 
simulation of the full process results in significantly different distributions of tbH='= selection variables, 
mainly in the distribution of the b-quark produced in association with the H^. 
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